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C/^ ' Abstract 

I propose a modification of the spherical infall model for the evolution of density 
fluctuations with initially Gaussian probability distribution and scale-free power spectra. 
I introduce a generalized form of the initial density distribution around an overdense 
I region and cut it off at half the inter-peak separation accounting in this way for the 

presence of neighbouring fluctuations. Contrary to the original predictions of the model, 
i resulting density profiles within virial radii no longer have power-law dependence on the 

On ' distance, but are well fitted by the universal formula of changing slope obtained as a 

. result of A^-body simulations. The profiles of halos of a given mass are in general flatter 

than the corresponding ones from the simulations, but the trend of steeper profiles for 
smaller masses is reproduced. The profiles of galaxy size objects are in general in better 
agreement with N-body simulations than those of larger ones. 
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^ • 1 Introduction 

■ 

High-resolution A^-body simulations with power-law initial power spectra suggest that den- 
^ ; sity profiles of dark halos in large range of masses are well fitted by a simple universal formula 

Pcrit,o (r/rs)(l + r/rs)2 
where Pcrit,o is the present critical density of the Universe and 5char is the characteristic density 

'^^'^^ = 3[ln(l + c)-c/(l + c)]' 

with V = 200. The scale radius Vs is related to the virial radius (the distance from the center 
of the halo within which the mean density is v times the critical density) by rg = r^/c and c is 
the concentration, the only fitting parameter in the formula. 

The density profile was observed to steepen from r~^ near the center of the halo to r^^ at 
large distances. This result seemed to contradict the prediction of the analytical spherical infall 
model (hereafter SIM) which for = 1 (the only case considered here) finds the profiles to 
be power-laws of the form 7--3{"+3)/(n+4) -^j-^gj^g fi jg ^j^g index of the initial power spectrum of 
density fluctuations. 



2 The modified spherical infall model 



I will argue here that the discrepancy between the two approaches is mainly due to the over- 
simplifications applied in the SIM. While such assumptions as spherical symmetry of the initial 
density distribution and the absence of peculiar velocities will be kept, the shape of the initial 
density distribution can in fact be made more realistic. This distribution is usually described 
by the expected overdensity within r[ provided there is a peak (overdense region) of height aa 
at Ti = 0, where a is the rms fiuctuation of the linear density field smoothed on scale R. If the 
initial probability distribution of fiuctuations is Gaussian and the filter is Gaussian the general 
form of this quantity as a function of Xi = r^/R can be found 
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{n + l)xf 




(3) 



This function is fiat near the center and only at large distances from the peak it approaches 
the Xj"''"'''^'* power-law applied in 

Although in the fiat Universe any overdense region bounds the mass up to infinite distance, 
in reality there are always neighbouring fiuctuations that also gather mass. As a way to emulate 
this conditions I propose a second modification of the initial density distribution in the form 
of a cut-off. One can think of two ways of estimating the cut-off scale. First, such scale could 
be found as a coherence scale of the overdense region defined by the expected overdensity (^) 
being equal to its rms fiuctuation. It turns out however, that a more stringent constraint is 
induced by the presence of other peaks (see 0). Therefore here the cut-off will be introduced 
at the half inter-peak separation Xi^p-p/2 for the most reasonable height of the peak, a = 3. In 
the case of n = — 1 we have Xi^pp/2=6.45. The generalized initial density distribution with a 
cut-off will be modelled by 

A. (^.\ - (^i(^O) (A\ 

^i,cnt\X^) ]^ _^ g(xi-Xi,pp/2)/«; 

with the width of the filter w = 1. 

According to the SIM the subsequent shells numbered by the coordinate x^ will slow down 
due to the gravitational attraction of the peak, stop at the maximum radius and then collapse 
by some factor / to end up at the final radius 

X-J[Ai^cnt{Xi) + 1] 

^— X 7 — \ • ("^J 

The simplest versions of the SIM adopt the value / = 1/2 motivated by the virial theorem 
and it will also be assumed here but a more realistic description can be found in p|. The final 
profile of the virialized halo is then 

^ -q + aa^)(l + ^0'f-V? (6) 



Pcrit,0 \XJ dx 

where g = ^r{t)/(j'^ is the correlation coefficient 



3 Comparison with the universal profile 

Since the measurements of halo properties from A^-body simulations were done at the state 
corresponding to the present epoch, the same condition will be applied for SIM calculations. 
Once the initial redshift zi is specified, equating the collapse time to the present age of the 
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Figure 1: The density profiles of tlie lialo of mass of order 3 x lO^^/i~^M0 for n = —1. The 
sohd hne shows the prediction of the SIM. The long-dashed one (NFWl) gives the result of 
the iV-body simulations with their fitted concentration, while the short-dashed curve (NFW2) 
presents the formula (P with concentration fitted to SIM results. 

Universe determines the overdensity of the presently virializing shell which ends up at the 
virial radius of the halo. When we adopt the normalization of the initial power spectrum 
(cTg = 1) and the conditions a = 3 and aa = 0.1 (for the linear theory to be valid) choosing the 
initial redshift Zi for a given spectral index n gives the comoving smoothing scale R with which 
the overdense regions are identified. The mass of the halo within the virial radius can then 
also be determined 

M=^Pcrit,o 7^ . (7) 



3 ' '"Vl + Zi. 

Figure |I| shows the density profile of a dark matter halo of galactic mass. The solid line 
presents the prediction of the SIM obtained from formula (j^) for n = ~1, zi = 600 and 
R = 0.188/i~^ Mpc. The final (virial) proper radius of the halo is = 0.231/i~^ Mpc and 
the mass M = 2.88 x IO^'^H'^Mq which correspond to a galactic halo. We see that the result 
of the SIM can be well fitted by formula (|l|) but the SIM profile is significantly flatter than 
the corresponding one from the simulations. The concentration parameters which measure 
the steepness of the profile (the higher c the steeper the profile) are c = 57.1 and c = 19.1 
respectively from the simulations and from the SIM. 

One of the main results of A^-body simulations was the dependence of the shape of the 
density profiles of halos on their mass. On the other hand, the standard prediction of the 
SIM gives the same profile independently of mass. However, with the improvements introduced 
above it is possible to reproduce the dependence of the profiles on mass. 

It is sometimes argued that if the density field is smoothed with a given scale R lower peaks 
end up as galaxies and higher ones as clusters. This, however, would violate the hierarchical 
way of structure formation since higher peaks collapse earlier. Another argument against such 
assumption comes from the calculations based on the improved SIM: the reasonable range of 
peak heights a between 2 and 4, which are most likely to produce halos, leads for a given 
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Figure 2: The dependence of the concentration parameter c on mass (given in units of the 
present nonhnear mass M^,). The dashed curves show the results of A^-body simulations and 
the solid lines represent the results of the SIM. In each set the curves from top to bottom 
correspond to n = 0, —0.5, —1 and —1.5 respectively. 

smoothing scale to the range of masses spanning only one order of magnitude, while in A^-body 
simulations halos with masses spanning few orders of magnitude are observed. This suggests 
that the dependence on mass should rather be related to the initial smoothing scale. 

The dependence of the shape of the profiles on mass obtained with these assumptions is 
shown in Figure 0. The solid lines give the values of the concentration parameter c obtained 
by fitting the formula (|l|) to the results of SIM for different power spectra and the dashed lines 
are the corresponding values from the simulations 0. The overall trend of steeper profiles for 
smaller masses is reproduced and the agreement between the two approaches is significantly 
better for smaller masses. 

The spherical infall model provides simple understanding of the dependence of the shape of 
the halo on its mass: smaller halos start forming earlier and by the present epoch their virial 
radii reach the cut-off scale that accounts for the presence of the neighbouring fluctuations; 
more massive halos form later and their virial radii are not affected by the cut-off scale, their 
virialized regions contain only the material that initially was quite close to the peak identified 
with the smoothing scale corresponding to the mass. 
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